function newVector = shipSimple(currentVector, shipPos)

maxVel = 10;
curVel = norm(currentVector);                       %current velocity

if(curVel == 0)
    newVel = rand(1)*2;                             %random velocity between [0,2];
    newDir = rand(1,2)*2-1;                         %completely random direction
    newVec = newDir/norm(newDir);
else
    newVel = -1;
    while(newVel<=0)
        newVel = curVel + (rand(1)*2-1);            %current velocity + random number in interval [-1,1]
    end
    normVector = currentVector/curVel;              %normalized vector
                                                    %current direction = inv tangent(y/x) in degrees
    if(normVector(1)==0)
        if(normVector(2)>0)
            curDir = 90;
        else
            curDir = 270;
        end
    elseif(normVector(1)>0)
        curDir = atand(normVector(2)/normVector(1));
    else
        curDir = atand(normVector(2)/normVector(1))+180;
    end
    
    newDir = curDir + (rand(1)*90-45);                 %current direction  + random number in interval [-45,45]
    if(newDir>=360)
        newDir = newDir-360;
    elseif(newDir<=-360)
        newDir = newDir+360;
    end
    newVec = [cosd(newDir), sind(newDir)];             %calculate the new unit vector
end

if(newVel > maxVel)
    newVel = maxVel;
end
% 
% newDir
% newVel

newVector = newVel*newVec;